Node-like excitations in superconducting PbMoeSg probed by scanning tunneling 

spectroscopy 
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We present the first scanning tunnehng spectroscopy study on the Chevrel phase PbMoeSg, an 
extreme type II superconductor with a coherence length only slightly larger than in high- Tc cuprates. 
Tunneling spectra measured on atomically flat terraces are spatially homogeneous and show well- 
defined coherence peaks. The low-energy spectral weight, the zero bias conductance and the temper- 
ature dependence of the gap are incompatible with a conventional isotropic s-wave interpretation, 
revealing the presence of low-energy excitations in the superconducting state. We show that our 
data are consistent with the presence of nodes in the superconducting gap. 
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The Chevrel phase family of superconductors attracted 
great attention in the 1970s and '80s since several of its 
members possessed unprecedentedly high upper critical 
fields while others displayed astonishing properties re- 
sulting from the coexistence of superconductivity and 
magnetismii Interest in these compounds faded some- 
what with the advent of high temperature supercon- 
ductors (HTS). However, in the quest to understand 
the pairing mechanism in HTS, research has recently 
diversified into other strongly correlated electron sys- 
tems which exhibit unusual superconductivity. For in- 
stance, heavy fermions^^i^ borocarbides^ii^ organics^ 
and ruthenates^i^ have already provided vital insight into 
the fundamental questions of gap symmetry and the in- 
terplay between superconductivity and magnetism. In 
this context, the properties of the Chevrel phases render 
them particularly relevant and interesting. Among these 
compounds, PbMogSs is one of the most remarkable with 
its high critical temperature of 14 K and a strikingly 
high upper critical field of about 60 T, corresponding 
to a very short coherence length ^ w 23 A (estimated 
from the critical fieldiS using Ginzburg-Landau theory). 
This value is comparable to those found in the HTS 
(e.g. 11, 20, 30 A for BisSrsCaCusOg+a, YBaaCgOT-a, 
La2_a;Sra;Cu04 respectively) With a k = A/^ above 
100, PbMoeSg is clearly an extreme type II supercon- 
ductor. 

Scanning Tunneling Spectroscopy (STS) today plays a 
central role in the study of novel superconductors at an 
atomic scale. In particular, STS investigations have re- 
vealed significant differences between classical BCS and 
high temperature superconductors notably the very 
large gaps in the spectra consistent with d-wave symme- 
try and the spatial dependence of the electronic structure 
in the vortex cores which differ radically from what is ob- 
served in conventional superconductors . One possible 
contribution to these differences is the very small size of 
the Cooper pairs in HTS, resulting in a deviation from 



mean field theory. Compared with HTS, PbMogSg is sto- 
ichiometric, its structure 3D and, unlike other Chevrel 
phases, is not known to display competition between su- 
perconductivity and magnetism. We consider it to be 
a bridge between classical typc-II superconductors and 
HTS. When the Chevrel phases were first studied, it was 
assumed during most investigations that they were clas- 
sical s-wave superconductorsi similar to other type II 
superconductors known at that time. However, the very 
small Cooper pair size raises questions over the micro- 
scopic nature of the superconducting state in PbMogSg 
and, more generally, the role of short coherence length in 
unconventional superconductivity. 

In this paper we present the first STS study of 
PbMogSs single crystals. We show tunneling spectra 
measured on atomically flat cleaved surfaces at different 
temperatures. The central result reported here is that 
the spectra obtained reproducibly on different crystals 
are clearly incompatible with the classical s-wave model 
and reveal the presence of low-energy excitations between 
the superconducting coherence peaks. 

The high-quality PbMogSg single crystals were pre- 
pared by both chemical vapor transport in silica tubes 
and high temperature melt techniques in sealed Mo cru- 
ciblesJ^ They display a sharp superconducting transition 
at Tc = 14 K with width AT^ = 0.3 K determined by 
ac susceptibility measurements. PbMogSg crystallizes in 
a "quasi-cubic" rhombohedral lattice, space group i?3,— 
with the three- fold axis along the [111] direction. We 
obtained flat (001) surfaces by cleaving the crystals us- 
ing a blade under an optical microscope. The STS mea- 
surements were carried out using a home-built scanning 
tunneling microscope, featuring an XY-stage^^ to target 
the micrometer-sized crystals. We use electrochemically 
etched Ir tips and measure the differential conductivity 
using a standard lock-in technique. 

The surface topography imaged at 1.8 K with the tip 
perpendicular to the (001) plane reveals large flat ter- 
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FIG. 1: Topographic characterization, (a) Large-scale topog- 
raphy (T = 1.8 K, Rt = 50 MQ). (b) Corresponding height 
histogram showing a unit step of height 5.9 ± 1.1 A. 



races separated by long parallel steps. Fig. [T] shows a 
typical step of height 5.9 ± 1.1 A. This corresponds to 
the distance between the Pb planes which are the natural 
cleavage planes of PbMoeSg. The surfaces investigated 
were atomically flat with a RMS roughness under 0.7 A. 

Homogeneous spectra have been obtained at a large 
number of positions on different crystals along lines 
(Fig. [2^) and on maps extending over a range of several 
hundred Angstroms. We have verified that the spectra 
obtained by varying the tunnel resistance Rt all collapse 
onto a single curve (Fig. [2)3) , which confirms true vacuum 
tunneling conditions. It is also important to mention that 
our measurements have always shown the same spectro- 
scopic signature regardless of sample and cleavage. The 
spectra present well-developed and homogeneous coher- 
ence peaks. The zero bias conductance (ZBC) is homo- 
geneous with a value of ^ 30% of the high-energy con- 
ductance (Fig. This unusually large value, also re- 
ported by early planar-junction experiments,— could in 
principle be due to an extrinsic influence such as surface 
impurities. However, the high reproducibility of the ZBC 
leads us to believe it to be intrinsic to this material. Fur- 
thermore, we stress that PbMogSs docs not suffer from 
surface contamination or ageing effects. 

The shape of the spectra reveals the presence of low- 
energy excitations which cannot be explained by thermal 
broadening and is reminiscent of that observed in HTS 
(Fig. [3]). We have therefore performed fits to our ex- 
perimental spectra using a single-band model and dif- 
ferent singlet-gap symmetry scenarii. STS is a direct 
probe of the local quasiparticlc DOS and is therefore not 
sensitive to the phase of the order parameter. Hence it 
cannot determine whether a sign change occurs, but it 
can reveal the presence of regions where the gap func- 
tion is zero (nodes) or nearly zero (node- like) . Using the 
BCS Green's function and including lifetime effects due 
to impurity scattering, the quasiparticlc dcnsity-of-states 
(DOS) reads: 



N(uj) = - i Vim 
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FIG. 2: (a) Spectroscopic trace along a 1100 A path with one 
spectrum every 8.8 A. The spectra are offset for clarity (T = 
1.8 K, Rt — 25 MQ). (b) Normalized conductance spectra 
obtained at various Rt from 20 to 50 MQ. (c) Histogram of 
the zero bias conductance in ^ 7500 spectra normalized to 
the conductance at 8 mV. The standard deviation is 0.04. 



scattering rate. We have calculated the DOS assuming a 
spherically symmetric Fermi surface; the relevance of the 
real band structure will be discussed below. 

The gap representations compatible with the rhombo- 
hedral symmetry are very similar to those for the hexag- 
onal group in Ref. [l^. We shall focus on the lowest-order 
gaps (s and d) which are listed in Table HI Representa- 
tions ^2 and ^4 give identical DOS as they only differ 
by a phase factor. Furthermore, representation J^^, al- 
though behaving differently close to the nodes, gives a 
DOS which is indistinguishable from that of .^2 and ^4 
once integrated over the Fermi surface. We thus group 
these three representations under the collective label "d- 
wave".— In contrast, the typical DOS from the last d- 



TABLE I: Lowest-order representations compatible with 
rhombohedral symmetry for the gap Afe = Aoij{k). z is 
chosen along the 3-fold axis. 



where = ^ ep are the band energies relative to 
the Fermi level, the superconducting gap and F the 
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FIG. 3: (a) Experimental dl /dV spectrum at 1.8 K (thick 
line) compared to the s-wave and d-wave fits described in the 
text (dashed and solid lines respectively), (b) dl /dV spectra 
obtained at various temperatures (thick lines) and d-wave fits 
(thin lines). Each experimental spectrum is normalized so as 
to ensure state conservation within the measurement energy 
range and offset by 0.45 with respect to the previous one for 
clarity. 



wave representation, ^5, shows greatly reduced peaks 
(similar to higher-order 5-wave cases not shown in Ta- 
ble HI which are totally incompatible with the observed 
spectra and will therefore not be further discussed here. 

The average experimental spectrum at 1.8 K is shown 
in Fig. [3^ together with the fits. The calculated spec- 
tra take into account temperature smearing and a Gaus- 



sian broadening with a half-width at half-maximum of 
0.4 mV. This value corresponds to the lock-in amplitude 
and sets the maximum experimental resolution. The fits 
are defined to optimally reproduce the peak energy posi- 
tions and widths. The isotropic s-wave for Aq = 2.6 mcV 
and r = 0.2 meV clearly fails to describe the spectrum 
between the coherence peaks. In direct contrast, the d- 
wave cases for Aq = 2.95 meV and F = mcV provide a 
much better fit in the low-energy region which is the most 
sensitive to changes in the gap symmetry. Attempting to 
fit the ZBC with a s-wave symmetry by raising F or in- 
troducing a normal residual DOS^ fails either due to a 
dramatic coherence peak reduction or a poor description 
of the spectrum at low energies. Furthermore, a mixed- 
symmetry gap function (e.g. 3£i -\- ^4) always results 
in a deterioration of the fit in the low-energy region. We 
emphasize that within the rf-wave scenario, the measured 
ZBC is a natural consequence of the low-lying excitations 
and thermal smearing. Note that the experimental coher- 
ence peaks are larger than those from d-wave fits. One 
possible explanation for this could be a spectral weight 
transfer from high energies to the peaks due to strong 
coupling effccts i^^i^^ 

Fig. shows spectra taken at different temperatures 
with their corresponding d-wave fits. We first observe 
that superconductivity measured at the surface is sup- 
pressed at exactly the bulk Tc = 14 K. This confirms that 
our measurements probe the bulk properties of PbMogSs. 
The d-wave models fit well at all temperatures, in con- 
trast to the s-wavc which consistently fails to reproduce 
the spectral shape at low energy. The gap amplitudes 
Ao(T) from the fits to our experimental spectra are plot- 
ted as a function of temperature in Fig. [51 together with 
the solutions of the BCS gap equation for both the s- 
wave and d-wave symmetries. These were calculated by 
selecting the simplest separable coupling interaction pre- 
serving the gap symmetry, i.e. Vk,k' = Vb'0(A;)V'*(fe')j 
and determining Vq so as to fix the zero-tcmperature gap 
at our experimental value (2.95 meV). The s-wave solu- 
tion gives a Tc of 19.4 K, as expected from the well-known 
2Ao/A:b?c ratio of 3.52 but clearly at variance with our 
experiment. On the other hand, the Tc from the 3D d- 
wave model is 13.6 K, giving a ratio of 5.0 in agreement 
with our experimental value of 4.9 ± 0.3 (Fig. [4^). These 
results are also shown renormalized by their respective Tc 
values in Fig. [Us, highlighting the facts that the presence 
of nodes affects the shape of the curve (and not only the 
2Ao/A:bT'c ratio) and the experimental gap closes faster 
with increasing temperature than the theoretical curves. 

The large low-temperature gap is also consistent with 
the high upper critical field: using the Fermi velocity 
from Ref. [23, up ~ 4 x 10"* ms^^, we obtain a rough 
estimate of the coherence length £ ~ w 28 A, in 
line with the value derived earlier from Hc2'^ This very 
small coherence length and the large coupling ratio of 4.9 
further confirm the similarity with the high- Tc materials, 
which has already been emphasized by.Ucmura et al.— 

Let us now turn to the effect of the true band structure 
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FIG. 4: (a) Temperature dependence of the superconducting 
gap from the spectra in Fig. ^jp (circles) compared to the 
s and c?-wave cases from solving the BCS gap equation for 
Ao(0) = 2.95 meV (dashed and solid curves respectively), 
(b) Same data with gaps rescaled by Ao(0) and temperatures 
rescaled by the respective Tc values. 

of this material. Our FLAPW— calculations, in agree- 
ment with previous calculationsfi reveal that two Mo-d 
bands (of width 0.7 cV) cross the Fermi level, both con- 
tributing roughly equally to the total DOS at ep. This 
justifies a single-band approach in calculating the super- 
conducting properties. The normal state DOS appears 
constant at the energy scale of our measurements. How- 
ever, the calculated Fermi surface (FS) is not spherically 
symmetric and exhibits significant nesting at Q ~ 0.62 
and 1.08 (27r/a) along the [111] direction. Most strik- 
ingly, one of the two FS sheets displays two lobes above 
and below the basal plane. This could be of importance 
in the case of highly anisotropic superconducting gaps. 
In effect, the presence of two bands and the FS topol- 
ogy influence the gap distribution and could thus explain 



the slight disagreement for the peak heights and ZBC 
between the c?-wave models and the observed spectra as 
well as the deviations from the calculated temperature 
dependence. The details of the Fermi surface could also 
influence the spectra through the tunneling matrix el- 
ement. However, since we tunnel from a direction 54° 
off the main (3-fold) symmetry axis, the STM tip effec- 
tively probes all relevant parts of the Fermi surface. Any 
further development of the model presented here would 
demand the knowledge of the radial dependence of the 
gap function, itself requiring a complete description of 
the pairing interaction. 

In summary, we have measured the low-temperature 
characteristics of the quasiparticle DOS of superconduct- 
ing PbMogSs for the first time using scanning tunneling 
spectroscopy. In light of these results, we conclude that a 
highly anisotropic gap function is the likely explanation 
for the low-energy excitations observed in our spectra. 
In particular, we have shown that our measurements are 
well described by 3D d-wave models compatible with the 
rhombohedral symmetry of PbMogSg- Such d-wave su- 
perconductivity in PbMogSs could result from the very 
short coherence length which would favor the appear- 
ance of a repulsive component in the coupling interaction. 
The combination of possible d-wave symmetry, a large 
2Ao/fcBT'c ratio and short coherence length in PbMogSs 
is strongly reminiscent of the HTS and calls for further 
investigation. We believe that vortex core spectroscopy, 
which has already revealed important differences between 
classical superconductors and HTS, would be the ideal 
tool to confirm the role played by an intrinsic short co- 
herence length in unconventional superconductivity. 
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